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We analyze the dynamics of a Bose-Einstein condensate undergoing a continuous dispersive imag- 
ing by using a Lindblad operator formalism. Continuous strong measurements drive the condensate 
out of the coherent state description assumed within the Gross-Pitaevskii mean-field approach. Con- 
■») | tinuous weak measurements allow instead to replace, for timescales short enough, the exact problem 

■ with its mean-field approximation through a stochastic analogue of the Gross-Pitaevskii equation. 

' The latter is used to show the unwinding of a dark soliton undergoing a continuous imaging. 



o 

Q 



> 
in 



o 

c3 



03.75.Fi,05.30.Jp,05.40.-a 



I. INTRODUCTION 



The interplay between quantum and classical descriptions of the physical world and the role of the measurement 
process are still at the heart of the understanding of quantum mechanics jjj . The related theoretical debate has been 
greatly enriched in the last decades by the realization of new experimental techniques aimed at producing quantum 
£C) , states without classical analogue, such as entangled or squeezed states, or to explore phenomena which are intrinsically 
C^l ■ quantum mechanical, like quantum jumps. All this occurred also having in mind practical implications, such as the 
improvement of sensitivity of various devices operating at or near the quantum limit |^,^| . 

Recently, the production of a novel state of matter - Bose-Einstein condensates of dilute atomic gases - has opened 
a new road to explore macroscopic quantum phenomena with the precision characteristic of atomic physics [Q . Bose- 
Einstein condensates are naturally produced by cooling down atomic gases at ultralow temperature with the phase 
transition occurring in the 100 nK - 1/iK temperature range, i.e. when the thermal de Broglie wavelength becomes 
comparable to the average spacing between the atoms of a dilute (peak density 10 13 — 10 15 atoms/cm 3 ) trapped gas. 
""^ ■ Usually samples of Bose-Einstein condensates are made by 10 3 — 10 s atoms in the case of alkali species and 10 9 
or more in the case of hydrogen j5j . Their intrinsically small heat capacity does not allow for a direct manipulation 
and probing with material samples, such as microtips or nanostructures, since the thermal contact with the latter 
will induce a sudden vaporization of the Bose sample. Thus, manipulation and probing of Bose-Einstein condensates 
has been achieved so far only by using light beams. These probes can be classified according to the resonant or 
nonresonant nature of their interaction with the measured atomic sample. In the former case, the condensate interacts 
with a laser beam resonant (or close to resonance) with a selected atomic transition. The output beam is attenuated 
proportionally to the optical thickness of the condensate (also called column density), i.e. the condensate density 
integrated along the line of sight of the impinging beam. The absorption of the photons leads to recoil of the atoms 
then strongly perturbing the condensate. For typical values of intensity and duration of the probe light the condensate 
is strongly heated and a new replica has to be produced to further study its dynamics. From the viewpoint of quantum 
measurement theory this measurement is of type-H since it destroys the state of the observed system and forbids the 
study of the dynamics of a single quantum system An alternative technique, called dispersive imaging, allows for 
repeated measurements on a Bosc-Einstcin condensate and, as an extreme case, its continuous monitoring. In this 
measurement scheme off-resonance light is scattered by the condensate, thereby locally inducing optical phase-shifts 
which can be converted into light intensity modulations by homodyne or heterodyne techniques, for instance by using 
phase-contrast j?J or interference techniques |8| . Since the laser beam is off- resonant the absorption rate is small and 
the heating of the condensate is accordingly small. Thus, multiple shots of the same condensate can be taken - a 
type-I measurement in the quantum measurement theory language - allowing the study of the dynamics on the same 
sample. This has allowed to overcome the unavoidable shot-to-shot fluctuations always present in the production 
of different samples of Bose-Einstein condensates. Several phenomena whose observation rely on imaging at high 
accuracy the condensate, such as its formation in non-adiabatic conditions [l^jhort (To| and long |ll[] wavelength 
collective excitations, superfluid dynamics Jl^ |, vortices formation and decay 
this technique. 
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The repeated nondestructive monitoring lends itself to another question: is the measurement process influencing the 
dynamics of the condensate? Answers to this question are important also for disentangling the intrinsic condensate 
dynamics from the artifacts induced by the underlying measurement process. As we will see in this paper, the effect 
of the measurement can also be intentionally amplified to allow for unusual manipulation of the condensate itself. 
Dispersive imaging, in its more idealized form, is preserving the number of atoms and therefore represent a particular 
example of quantum nondcmolition (QND) measurement |^,|l^,^6|. We know that even quantum nondemolition 
measurements affect the state of the observed system, their unique feature being that the nondcmolitive observable 
maintains the same average value, albeit the probability distribution of its outcome can be affected as well as the 
average values of all the conjugate observables. Thus, we do expect that during the nondemolitive measurement of 
the condensate atom number there will be a measurement-induced nontrivial dynamics for both the variance of the 
monitored observable and the average values of all conjugate observables. Besides gaining insight into the dynamics 
of the measurement process, our model allows for the description of the irreversible driving of the condensate toward 
nonclassical states. 

In this paper we try to answer the question formulated above by building a realistic model for dispersive mea- 
surements of the condensate, extending the results reported in fl^l to the weak measurement regime. The plan of 
the paper is as follows. In Section II we introduce the dispersive coupling between atoms and light and derive the 
reduced master equation for the Bose condensate tracing out the variables of the electromagnetic field degrees of 
freedom. Under controllable approximations we obtain a Lindblad equation which allows to estabilish the rates for 
phase diffusion and depletion of the condensate during the dispersive imaging. In Section III we unravel the Lindblad 
equation in terms of a stochastic differential equation that, in the unmeasured case, corresponds to the description 
of a single N-body wavefunction. A solution of the stochastic N-body equation is discussed in the limit of strong 
continuous measurement, leading to the squeezing of number fluctuations, the main result described in |T^ ]. In the 
opposite limit of weak measurement and for an initial mean-field state this stochastic equation becomes the stochastic 
counterpart of the Gross-Pitaevskii equation, as discussed in Section IV. Its limit of validity is discussed in Section 
V by comparing its evolution for various parameters versus the exact evolution in the simple situation of a two mode 
system schematizing a condensate in a double well potential. This allows us to analyze, in Section VI, the effect of 
the measurement on the evolution of a condensate initially prepared in a soliton state. More general considerations 
on the potentiality of such an approach and its consequences are finally outlined in the conclusions. 



II. MASTER EQUATION FOR DISPERSIVE IMAGING OF A BOSE CONDENSATE 



Our main goal in this Section is to include the atom-photon interaction present in the dispersive imaging of a 
Bose-Einstein condensate into its intrinsic dynamics. First attempts in this direction have been discussed in the 



prototypical situation of a two mode condensate in p8 19 1. Let us start the analysis with the effective interaction 



Hamiltonian between the off-resonant photons and the atoms, written as 

Hint = ^ j ' d?x n(x) : E 2 :, (1) 

where n is the density operator of the atomic vapor, and E the electric field due to the intensity / of the incoming 
light. The coefficient xo represents the effective electric susceptibility of the atoms defined as xo = A 3 <5 /2tt 2 (1 + S 2 ), 
where we have introduced the light wavelength A and the light detuning measured in half-lincwidths T/2 of the atomic 
transition, S = (u — w at )/(F/2). 

We express the electric field in terms of creation and annihilation operators. In the Coulomb gauge it takes the 
form 



Eix. I'<=i^2\l ^^(a k e~ lckt+lWx - ale lckt ~ tkx ), (2) 



where Wk = c]k|, [ak, o^i] = ^k,k', and L 3 is the quantization volume. 



Equation (Of) allows us to write the reduced master equation for the atomic degrees of freedom by a standard 
technique, i.e. by tracing out the photon degrees of freedom ^0|. The decoupling between the two relevant timescales 
for the photons (settled by the the lifetime of spontaneous emission, of order of tens ns) and for the atoms (related 
to the oscillation period in the trapping potential, of the order of ms), allows to use the Born-Markov approximation. 
Thus we get the master equation for the reduced density matrix p of the condensate that, in the interaction picture, 
is written as 
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% = JtMpW ® Pr(*), flint] - ilhtlfftotW, [ / dt'H- m t(t'), p(t) ® p R (t)]]. (3) 

The last term on the righthandside contains two different contributions both of Lindblad type, Lip and L%p. As 
we will see soon, the former preserves the number of atoms in the condensate and is responsible for phase diffusion 
phenomena, while the latter changes the number of atoms leading to its depletion. 

Let us first concentrate on phase-diffusion which is a number conserving mechanism. To calculate it we insert the 
interaction Hamiltonian into the last term of Eq.(^), with n the condensate density operator. By introducing the 
Fourier transform of the density operator such that 



-( x ) = £ eiqx 7=' ( 4 ) 



we obtain 



n(k' - p')rc(k - p)p{a^apfa^ap) - n(k' - p')ph(k - p)(aj c a p aj c ,a p /) 
-n(k - p)pn(k' - p')(a' k ,a p 'a' k a p ) + pn(k - p)n(k' - p')(aj t .a p aj c ,ap') , (5) 

where (...) = Ttr[/Or . . .]. The photons are assumed to be in a coherent plane wave state with momentum along the 
impinging direction, corresponding to a wavevector k = k$z orthogonal to the imaging plane x — y. Hence 

L 3 I 

(al,a p ) = -^-5 k ,, P 5 pM , (6) 

where we have written the mean numbers of photons in mode ko in terms of the intensity of the incoming beam. 
Expressing the expectation values in normal ordering and using the fact that phase-diffusion processes conserve the 
number of bosons in the condensate, it is possible to show that all normal ordered expectation values involving four 
operators cancel exactly, obtaining 

LlP = ^mT (^)V d3fc<s d k + k °l- fc °) [»(-k),[ft(k),p]], (7) 

where we have used the continuum limit ^ k — > (L/2ir) 3 J d 3 k. Unless tomographic techniques are used as for 
instance in |2l]| , the image results from a projection of the condensate onto the x — y plane, by integrating along the 
z direction. This demands to project the dynamics of the condensate into the imaging plane. In order to write a 
closed 2D master equation to describe the x — y dynamics we assume the condensate wavefunction to be factorizable 
as tp(x,y,z) — (f>(x,y)A(z). Such factorization holds if the confinement in the z-direction is strong enough to make 
the corresponding mean-field energy negligible with respect to the energy quanta of the confinement, i.e. hio z » gp, 
where g = Airf^a/m, with a the s-wave scattering length, p the condensate density, and uj z the angular frequency of 
the confinement harmonic potential along the z direction, as recently experimentally demonstrated in p2| . We write 
n(k) = h(kj_)n(k z ) and we will use a Gaussian ansatz for the density profile along the z direction, namely 

I n f 2 fc 2 

n(k z ) = ] j 1 -e~^, (8) 

where £ is the lengthscale of the condensate in the z direction, the width of the Gaussian state A(z) under the 
abovementioned approximation. The effective 2D nonlinear coupling strength is g2D — 9 J dz\A(z)\ = g\pHj^. 
Consequently 



Lip = (27) / d2k± [ " ( ~ k±) ' [ " (k±) ' p]] (^) / dkz " 2 ( fc ^ti k + k °i - k o}- 



(9) 



By assuming that the typical length of the BEC in the z direction is much larger than the wavelength of the incoming 
laser, i.e., £ 3> A, we can calculate the value of the last integral, and it is equal to exp(— £ 2 k^/4fco). Our final result 
for the phase-diffusion contribution to the reduced master equation in the imaging plane is written as 
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lap = J d 2 ri J d 2 r 2 K(ri - r 2 ) [n(r x ), [n(r a ), p]], (10) 
where n(r) = ^^(r) is the 2D density operator, and K is the measurement kernel 

K (r) = J cP k exp(-£ 2 fc 4 /4fc 2 + ikr). (11) 

Eq. ( |l0| ) preserves the total number of atoms, and corresponds to a quantum nondemolition coupling between the 



atom and the optical fields 18 19|.|23| g5[. If the measurement kernel were a local one, K (y\ — r 2 ) ~ 5{v\ — r 2 ), Eq.(|l( 
would reduce to a Lindblad equation for the measurement of an infinite number of densities n(r). This assumes that 
no spatial correlation is established by the photon detection. However, the ultimate resolution limit in the imaging 
system depends on the photon wavelength, regardless of the pixel density of the detecting camera. The resolution 
lengthscale follows from Eq. ( |TT| ) as a width of the kernel 

Ar = (27^0 ^ = «A) 1/2 , (12) 

the geometrical average of the light wavelength and the condensate thickness £. Eq. ( |l0| ) can then be rewritten as 
L\p = Ji[n, [n, p]], where 71 is the phase diffusion rate, given by 

ll = ^L( d 2 ke^\m\ 2 \<t>(~n 2 - (13) 



(14) 



2hc 

We estimate its magnitude assuming a Gaussian profile in the x-y plane, i.e., |(/>(k)| 2 = e~ a k / 2 , obtaining 



71 = ~^r e £ 



1 - Erf 



a 2 k, 



2hca 2 ' 



where the last step holds for a well localized condensate, i.e. a 2 fco/£ 1. 

Let us now calculate the depletion contribution to the master equation for the condensate. We split the field 
annihilation operator into a term describing the condensate and another associated to the non-condensed particles, 
i.e., if) = ^c+^nc- We shall assume that the non-condensed particles belong to the continuum, so that their spectrum 
is the one of a free-particle hQq — Ti 2 q 2 /2m. Indeed, photons have large momenta with respect to the momenta of 
trapped atoms, so even if a small percentage of the photon momentum is absorbed by the atom, the atom is promoted 
into a high energy, unbounded state. Let us expand both field operators in terms of annihilation operators as 

tpc(x,t) = 6 c <Mx,£) 
^ NC (x, t) = ^=Y< e- ifM+iqx K (15) 

where cf>c is the condensate wave function. We get 
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Lip = -^-^ V^^p Y y w k w p^( n q' - ck> + c p') x 

kpq k'p'q' 

{e u ^~ ck+ ^Mq+ P -k)4>* c ( q ' +p'- fc')Tr R [at ak 6 q 6 c , [akvMc.P]] 

+e -im«-ck + c P )fe {q + p _ k) ^ c{q , +pl _ k > )Trn[a la p b^bt [a^a^b^p]]} . (16) 

Here the trace is taken over the reservoir of the condensate, which in this case consists of the non-condensed particles 
and the photons. The above expression contains depletion processes, in which a photon interacts with a particle in 
the condensate and, as a result, that particle is kicked out of the condensate. It also contains feeding processes, in 
which the reverse mechanism may take place. When one assumes that in the initial state of the non-condensate plane 
waves are empty, only the depletion process is relevant. In this hypothesis we find 



LiP = 72 + \ { 6 c 6 c, p}) , (17) 

where 72 is the depletion rate 
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72 = 5tK^) / ' d V 3 9^p l<Mq + p-k )| 2 5^-cko + cp). (18) 

We evaluate 72 in the thick condensate limit (£ 3> fc^ 1 ) by approximating |^c(q + p — ko)| 2 « (27r/L) 3 <5(q + p — ko), 
and then, by using the fact that Sl q <C cko, we obtain 

72 = a r ■ (19) 
snnc 

PromEqs.(0@ we see that the depletion rate is much bigger than the phase diffusion rate 72/71 = a 2 fco/47T 3 ^> 1, in 
accordance with Ref. [^5j . We estimate the magnitude of both rates using the following parameters for the condensate 
and its imaging, relevant for the case of 87 Rb: A = 780nm, \o — 10~ 23 m 3 , laser intensity / = 10~ 4 mW/cm 2 , and a 
typical size in the x-y plane of a — 10^m. Then 71 = 10 _6 s _1 and 72 = 10 s , corresponding to phase diffusion 
and depletion times of t\ = 7-f 1 = 10 5 s and t% = = 10 4 s, respectively. Although the depletion rate is larger that 
the phase diffusion rate, this last process can dominate because of their different scaling with the total number of 
condensed particles in the master equation. Indeed, the first is linear in AT, whereas the second one is quadratic, so 
for large number of particles (such as N — 10 7 ) phase diffusion occurs on a faster timescale than depletion. For this 
reason we will focus in the following on the phase diffusion contribution alone. 



III. STRONG MEASUREMENT: MEASUREMENT-INDUCED NUMBER SQUEEZING 

We will consider in the following the effect of strong measurements on the quantum state of the condensate. This 
has already been described in detail elsewhere [ p~7[ , thus here we only summarize the main results and their link to 
the following considerations. By neglecting the depletion term, the master equation for the condensate takes the form 
p = (— i/K)[H, p] — Lip. This equation preserves the total number of atoms in the condensate and corresponds to a 
quantum nondcmolition measurement of the atomic density via the optical fields. In order to get an insight into the 
equation, we introduce a two-dimensional lattice with lattice constant set by the kernel resolution Ar. In this way 
we get 

-Sj2[ni,[m,p}]. (20) 

Here "J; is an annihilation operator and n; = ^J^i is the number operator at a lattice site I. The frequency of hopping 
between any nearest neighbor sites (k,l) is uo ~ h/2mAr 2 , which is TtT 1 times the characteristic kinetic energy. The 
potential energy operator is V = J2i(^i n i + Gn 2 ), where Ui is the trapping potential and G = g2n/ Ar 2 . The effective 
measurement strength is S « J d 2 rK(r)/ Ar 2 — (2ir / Ar) 2 (irxokol /2hc) . In order to solve this equation we use an 
unraveling in terms of pure states such that p = \^f)(^\, the overline denoting the average over the unraveling 
stochastic realizations. The pointer states of Eq. (p0[ ) are not changed by the chosen unraveling ^6|. The pure states 
can be expanded in a Fock basis per site, \*S?) — zTfiV;} ^{JVjK-^i}) an d the amplitudes ^{Ni} satisfy the following 
stochastic Schrodinger equation (written in the Stratonovich convention) 

^{Ni} = ~tY1 h {Ni,N[}ij{N[} ~ ^V{Ni}^{N t } + ^{Ni} ^ [~ S ( N l ~ n lf + Sa l + W ~ n O l] > ( 21 ) 



where the homodyne noises have averages 9i(t) = and 9i t {t\)6i 2 (t 2 ) — 2S5i 1 j 2 8{t\ — 1<}). Here ni — X^at,} -^^{at,}! 2 * 
a 2 = J2{n,}(Ni — ni) 2 \ip{ Nl }\ 2 , h is the matrix element of the hopping Hamiltonian and Vjjv,} = J2i(UiNi + GNf) 
is the potential energy. When there is no hopping term (h = 0), an exact solution to this equation as a product of 
Gaussian-like wavefunctions is written as 

V^,W=e--e-^a t n^g^. (22) 

The population mean value per site ni(t) does a random walk, while its dispersion decreases as of (t) — of (0)/[l + 
4af(Q)Si\. Here of (0) is the initial dispersion in the number of atoms per site, and it scales as N for an initial 



dp 



1 

h 



-huj^2^l^ k + v,p 



(W) 
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coherent state. Thus, the measurement drives the quantum state of the condensate to a Fock state. When tunneling 
between different lattice sites is allowed (h ^ 0), localization in a Fock state is inhibited due to a competition between 
the measurement, which drives localization, and hopping, which tries to drive the state of the condensate towards 
coherent states. When measurement outweighs hopping, the final state of the BEC is a number squeezed state. The 
timescale in which squeezing is achieved is given by 

ts Q = (23) 

and the associated dispersion in the number of atoms per site for the asymptotic squeezed state is 

m = (conJS) 1 / 4 . (24) 

In a coherent state the number fluctuations per site are Poissonian, 07 = n~J 2 . The state is squeezed when (toni/S) 1 ^ < 
1/2 

n ; , i.e. sub-Poissonian atomic number fluctuations. This condition defines the strong measurement as 

mS 



> 1 . (25) 



Number squeezing of a Bose condensate has been experimentally observed in an optical lattice in |27[ . Our situation 
has an important difference from the latter case: since the squeezing is driven by the (Lindblad) measurement term 
the evolution into such states is irreversible, even after removal of the imaging photon field. From this viewpoint our 
squeezing technique is similar to the spin squeezing through quantum nondemolition measurements proposed in ]2§| ] 
and demonstrated in [^9). Of course, the system will eventually drift towards coherent, classical states due to the 
interaction with the external environment and the related decoherence ]3(| , for instance due to the thermal component 
or residual background gas in the trapping volume. 



IV. WEAK MEASUREMENT: STOCHASTIC GROSS-PITAEVSKII EQUATION 



Unraveling the Lindblad equation derived above leads, in the limit of weak measurements, to a stochastic equation 
for the condensate wavefunction. In the mean-field approximation this equation becomes the analogue of the (deter- 
ministic) Gross-Pitaevskii equation for the unmeasured system. From Eqs. (|To| , |TT| ) and ignoring the depletion term 
— L 2 p we obtain a continuum version of the master equation 



dp 
~dt 



h 



[H,p] 



d 2 n / d 2 r 2 K{r 1 -v 2 ) Kn),Kr 2 ),p]], 



(26) 



where H is the self-Hamiltonian of the system of N atoms in a 2D trap, 



H 



d 2 r 



2m 2 



(27) 



A nonlinear stochastic (Ito) unraveling of the master equation (|2q ) is 



d|¥) 



dfcff|¥) - dt / d 2 r x \ d 2 r 2 K{r 1 - r 2 ) An(r x ) An(r 2 ) |*) + / d 2 r dW{r) An(r) |* 



(28) 



Here An(r) — n(r) — (W|n(r)|W) and the Gaussian noises have correlators dW(r) = and dW{r\)dW{r2) — 2dtK(v\ — 
r 2 ). This unravelling corresponds to phase contrast measurement of the density of the condensate. The evolution of 
\^} given by Eq.(|2g|) describes a single realization of the experiment [^o| , ^ll -]33]| . 

For a single atom, N = 1, described by a wavefunction <f)(t,r) — (r|W), the stochastic Schrodinger equation takes 
the form 



d0(r) 



1 , 
--dt 
n 



rv2 

2m 



U(r) 



(r) + 



dW(r) 



d 2 r' |0(r')| 2 dW{v') 



(r) + C.T. 



(29) 



The second term in brackets follows from the last term in Eq.(p8|), while C.T. (counterterm) comes from the second 
term in J2S|). The latter is necessary to conserve the norm J d?r |^(r)| 2 = 1, and it is given by 



G 



C.T. 



dt 



K(0) + 2 J d 2 n K(r - n) |0( ri )| 2 - / d 2 n y"d 2 r 2 |^(r 1 )| 2 ^(r 1 -r 2 ) |0(r 2 )| 5 



(r). (30) 



This counterterm is badly nonlocal. Fortunately, to implement the stochastic terms numerically one can use 



d\cj)(r)\ 2 = 2dt 



dW{r) 



d 2 r' |0(r')r dW(r') 



dW(r) - / dV |<J>(r')| 2 dW(r') 



|^(r)| 2 + [^(r)(C.T.)+c.c] 



2 dt 



dW{v) 



d 2 r' \(j)(r')\ 2 dW(r') 



m\ 2 



(31) 



where the C.T. is used to cancel out the average of the stochastic terms squared. The stochastic terms affect the 
modulus of 4>{r) only, so Eq. (J3l|) is all that one needs to implement them. Eq.(|3l]) manifestly conserves the norm. 

For a weak measurement, when the squeezing of the quantum state is small, one can assume that the stochastic 
conditional state of N atoms in Eq. (E8h is a product mean-field state 



1 



*N\ 



d 2 r (f>(t, r)¥(r) 



N 



10), 



(32) 



with all the JV atoms in the condensate wavefunction <j>(t,r). The latter evolves according to an Ito stochastic 
Gross-Pitaevskii (SGP) equation 



d(f>(r) 



-dt 



2m 



V 2 + U(r) + (JV — 1) cj2D |0(r)| 2 



(r) 



dW{v) - / d 2 r' \4>{y')\ 2 dW(r') 



C.T. 



(33) 



In comparison to the case of JV = 1 (see Eq. (|29|) ) , this SGP equation has the usual extra nonlinear term (JV — 
l)52D|</>(r)| 2 but the stochastic terms containing dW are the same as for a single atom described by Eq.(|29|). This 
may seem strange because a state of JV 3> 1 atoms in the same quantum state can be found much faster than a state 
of a single atom. One might expect the stochastic terms, which describe backaction of the measurement on the state 
of the system, to become stronger with increasing JV. The backaction on the JV-atom state is indeed stronger than 
on the single atom state but its effect is divided over JV > 1 atoms instead of one. The backaction effect on the 
mean-field state </>(r) of each of the JV atoms is the same as the backaction on the quantum state of a single atom. 
Thus the stochastic terms for JV > 1 are the same as for JV = 1. 

As a final remark, we mention that a stochastic Gross-Pitaevskii equation has been proposed and studied for a 
quite different goal, namely to describe single trajectories through unraveling of the exact N-body quantum evolution 
of a Boson system |34[ ]. In the latter case the interpretation of the underlying stochasticity is obtained in terms of 
the randomness attributable to each quantum trajectory, to be confronted with the stochasticity that in our case is 
instead due to the opening of the condensate to a particular environment, namely the measurement apparatus. 



V. SGP VERSUS EXACT QUANTUM EVOLUTION: MEASUREMENTS IN A DOUBLE WELL 

In this Section we want to test the validity of the SGP equation in a significant but simple situation. To this 
end we consider the double well problem in the two mode approximation, and compare the quantum dynamics with 
measurement backaction with the dynamics given by the SGP equation. The Hamiltonian of the model is 

G 

H = e{a\ai + a\a 2 ) — hu!(a\a 2 + a\a\) + — [{a\) 2 a\ + [a^fa 2 ,], (34) 

where e is the mode frequency (assumed to be the same for both modes) , u> is the tunneling angular frequency, and G is 
the two-particle interaction strength. We perform phase-contrast imaging on each site, and we assume that the kernel 
resolution is much shorter than the distance between sites, so that there is no cross term due to the measurement. 
The Ito version of the stochastic Schrodinger equation for the state ket reads 
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d|* Q ) =dt 



■^-§(«i-<m)) 2 -§(n 2 -<™2>) 2 



dWi(m - (m))l*q) + dW 2 (n 2 - (n 2 ))\V<. 



(35) 



where the noise satisfies dW a = and dW a dWp — 2K a pdt, with a, (3 = 1,2, and lf aj 8 = S5 a ^ the measurement 
kernel. By assuming a total of iV atoms, distributed between the two sites, N = N\ + N 2 , we can expand the state 
in terms of Fock states at each well 



JV 



I* 



(36) 



fe=0 



We solve numerically the corresponding equation for the coefficients ipk(t) starting from a mean field state with 
equal mean populations in each well, |^)t=o = ~^j(^ a \ + ~^ a 2) N \ty- The equation provides us with the full 
quantum evolution including the measurement backaction. We want to compare it with the one that results from the 
stochastic GP equation. The GP state is 



|*gp> = 



1 



[fa(t)a\ 



Mt)4] N \o), 



(37) 



where the wavefunctions fa and fa satisfy the following set of Ito SGP equations, that follow from Eqs.(J30|,|33|) 



#i 



: dt[(e + (N-l)G\fa 



ufa] + dWi(l - \fa\ 2 ) + Sdt[-1 + 2\fa\ 2 - (\fa\ 4 + |0 2 | 4 )], 



h = + ( N ~ l)G\fa\ 2 )fa - ufa] + dW 2 (l - \fa\ 2 ) + Sdt[-1 + 2\fa\ 2 - (Ifal* + \fa\% 



(38) 



and we take the same initial state as in the quantum evolution, namely a coherent state fa(0) = fa (0) = l/v2 with 
balanced populations ni(Q) = n 2 (0). In the Fig.l we compare the time evolution of the quantum and SGP mean 
populations in one well for increasing values of the measurement strength. In these simulations the total number of 
particles was TV = 100, and the nonlinear coupling was G = Q. In this case, the Hamiltonian involves only one-body 
terms, so the mean field evolution (based on coherent states) must exactly coincide with the quantum one for the 
case of zero measurement (S = 0). For small measurement strengths, nS/to <C 1 (here n denotes the average number 
of particles per site), the quantum state of the condensate is still, to a high degree of accuracy, a coherent state, so 
the SGP evolution and the quantum one coincide. In Fig. la we see that the agreement is good even for times much 
larger than 1/nS. This timescale is relevant for the strong measurement case of the previous Section, since it sets the 
time after which an asymptotic number squeezed state is reached. As we increase the measurement strength and reach 
nS/cu > 1 (see Fig. lb), the mean number of particles per well given by the SGP and the quantum evolution depart 
appreciably, not surprisingly since such strong measurements squeeze the quantum state of the condensate driving it 
outside the description in terms of coherent states, i.e. the associated basis for the Gross-Pitaevskii equation. 




FIG. 1. Comparison between the SGP and the exact quantum evolution for weak measurements on a non-interacting Bose 
condensate. The population n(t) in a given well is plotted as a function of time for a single stochastic realization, expressed 
in units of cycles of the Rabi oscillations. The total number of particles is iV = 100. (a) The SGP evolution and the quantum 
one coincide for an effective measurement strength nS/u) = 10 _1 even for times much greater than 1/nS = 10/w. (b) The 
departure of the SGP dynamics from the exact one at the time 1/nS = 1/w is evident for a stronger measurement coupling 
nS/u = 1. 
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The Gross-Pitaevskii evolution can depart from the quantum one not only due to the measurement backaction but 
also due to the nonlinearity of the interactions. In Fig. 2a we show the SGP and quantum evolutions for S/u) = 1CP 3 
and G/uj — 10~ 3 , corresponding to the same initial state as in previous figures. We see that the inclusion of the 
nonlinearity causes the SGP evolution to depart from the quantum one. To gain further insight we introduce a 
quantity which measures the depleted fraction of atoms from the best mean field state, i.e. a mean field state that is 
the closest to the exact quantum state. Its definition is 

min / I \ 

D={A,<j>}[\--(* Q \Jc\* Q ) > ) (39) 

where the operator c is c = \/~Aa\ + e^^/l — A 2 a 2 - In Fig. 2b we plot this depletion for the same simulation of Fig. 
2a. The depletion in Fig. 2b is small, less than one atom is depleted from the condensate, but, as we see in Fig. 2a, the 
SGP evolution departs from the exact evolution. This departure is attributable to the inclusion of the nonzero G. 




FIG. 2. Comparison between the SGP and the exact quantum evolution for a weakly interacting {G/uj = 1CP 3 ) Bose conden- 
sate subjected to a weak measurement nS/ui = 10 _1 . (a) The population difference An(t) = 712 (t) —ni(t) is plotted for a single 
stochastic realization. Initially the condensate is in a coherent state with equal populations in each mode, ni(0) = ri2(0) = 50. 
(b) Depletion from a mean field state for the same realization as in (a) . 

In Fig. 3a we show the GP and quantum evolution for S = and G/u> = 10 -3 . Since no measurements are 
performed, an initial balanced population (ni(0) = ^(0)) would remain balanced for all times, both at the GP and 
quantum level. For this reason we take an initial unbalanced population, ni(0)/n 2 (0) = 3/2, which due to the hopping 
term triggers Rabi oscillations between the two wells (just as the measurement did when we took initial balanced 
populations). It follows from the figure that the GP dynamics departs from the quantum one. We plot in Fig. 3b 
the depletion corresponding to Fig. 3a. Again, this depletion remains small, less than one atom is depleted from the 
condensate. 




cot cot 
FIG. 3. Comparison between the standard GP and the exact quantum evolution for a weakly interacting (G/oj = 1CP' 5 ) 
unmeasured Bose condensate, (a) The population difference An(t) = 7i2(i) — n\{t) is plotted as a function of time. Initially 
the condensate is in a coherent state with unbalanced populations, ni(0) = 60, 712(0) = 40. (b) Depletion from a mean field 
state for the same realization as in (a). 

As we can see in Figs. 2, 3, both the SGP and GP evolutions depart from the exact evolution even in the weak 
measurement and interaction limits when the depletion from the condensate is small. The derivation of the backaction 
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terms in the SGP equation requires only one assumption: all the atoms are in the condensate. In contrast, the 
derivation of the interaction terms (both for SGP and GP equations) not only assumes that all the atoms are in the 
condensate but also makes further approximations to describe the evolution of the condensate wavefunction (p. This 
is why G ^ causes both the SGP and GP evolutions to depart from the exact quantum evolution even for negligible 
depletion. However, even after the departure the SGP and GP evolutions remain qualitatively similar to the exact 
evolution, see Figs. 2a and 3a. 

To summarize, the accuracy of the SGP equation is limited by the weak measurement condition, nS/u < 1, and by 
the nonlinear interaction. For a weak measurement the accuracy of the SGP equation is the same as that of the GP 
equation. In the next Section we apply the SGP equation to show how the measurement can trigger the unwinding 
of a dark soliton [E35|3(| in a realistic experimental setup. 



VI. MEASUREMENT-INDUCED UNWINDING OF A DARK SOLITON 



Once delimited the validity of the stochastic Gross-Pitaevskii equation in the simple situation of a two-mode system 
we can apply it to the more complex case of imaging of a condensate state with a nontrivial phase such as a soliton. 
Even in the limit of weak measurement, with the solution still approximable in terms of Gross-Pitaevskii coherent 
states, the effect of the measurement is present and affects the observable conjugate to the atom number, i.e. the 
phase of the condensate. As an example of application of the SGPE let us consider an isotropic harmonic 2D trapping 
potential V(r) = (rati 2 /2)(x 2 + y 2 ). The condensate is assumed to be in the Thomas- Fermi limit of strong repulsive 
interaction where the ground state wavefunction can be well approximated by 



^cs(r) = Jl^m • (40) 

V N 92D 



The constant /i = ClWNg^mjTr, the chemical potential, is chosen so that the wavefunction is normalized to 1. Let 
us use as the initial state a dark soliton |55| imprinted on the Thomas-Fermi ground state 

*-"■'> --"(f) ^ss^- < 41 ' 



Here I — 0.6/im is the healing length at the peak density in the ground state (|40|). 

In our numerical simulations we assume the following parameters relevant for 87 Rb: mass m — 1.4 x 10~ 25 kg, 
scattering length a = 5.8nm, xo = 10 _23 m 3 , wavelength A = 780nm. The width of the gaussian A(z) is assumed to 
be £ = 10/mi. With these parameters the resolution of the kernel is Ar = 5/.tm. We assume a 2D harmonic trap 
frequency £1 — 2ir x 10s _1 . With N = 5 x 10 5 we get a 2D Thomas-Fermi radius of 31/im. We also assume the laser 
intensity of I = 10~ 4 mW/cm 2 . 

The former parameters give a weak measurement strength in the sense discussed in Sections III and IV, so that the 
use of the SGP equation is justified. Indeed, this continuous problem can be mapped on the lattice model Eq.(^o|), 
where the lattice constant is the kernel resolution Ar. Given the Thomas-Fermi radius of 31/im, and Ar = 5/im, 
we estimate the lattice to be composed of 100 sites with an average of n = 5000 atoms per site. The measurement 
strength is equal to S — 7 x 10 _5 s , and the effective hopping frequency for the lattice model is u> = 14s -1 . Therefore 
nS/uj «2x 10~ 2 <C 1, thus confirming a weak measurement regime. 

To simulate the continuum SGP equation ( |33| ) we discretized it using a lattice constant 2tt times smaller than the 
kernel resolution Ar. The program uses a split step method: the fast Fourier transform was used to carry out the 
time integration of the kinetic term and of the nonlocal terms involving the kernel, and the potential and nonlinear 
coupling terms were integrated in time in the position representation. A cross-section along the x-axis through the 
probability density of the initial state ([D]) is shown in (a) of Fig. 4. Cross sections through probability densities at 
later times after the probe light beam has been sent on the condensate are shown in cases (b) and (c) of Fig. 4. For 
comparison, the time evolution without the measurement does not result in any soliton unwinding. 
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FIG. 4. Soliton unwinding during a continuous imaging. Figures show the cross section of the 2D probability density 
\§(x, y)\ 2 of the Bose condensate along the a;-axis at t = ms (a), t = 50 ms (b), and t = 200 ms (c) since the beginning of the 
continuous imaging. The dark soliton is greying progressively and in the third profile there is no trace of the initially imprinted 
pattern. The laser intensity is I = 10~ 4 mW/cm 2 , and the condensate wavefunction is normalized to unity in 2D. 



The soliton unwinds after roughly 50 ms. This time is much shorter than the depletion time £2 = 10 4 s discussed 
at the end of Section II. The measurement induces soliton unwinding much earlier that any detectable depletion of 
atoms occurs. Figure || suggests that the unwinding will manifest itself by filling up the soliton core with atoms. Such 
a greying of the dark soliton can also occur through a different mechanism that involves collisions between condensate 
atoms. In Ref. |3lj] it was demonstrated that the dark soliton can grey on a timescale of tens of ms because its core 
fills up with non-condensed atoms (quantum) depleted from the condensate as a result of atomic collisions between 
condensed atoms. This result is supported by Ref. p§| where it was shown that the quantum state with minimal 
depletion has depletion strongly concentrated in the soliton core. This quantum depletion process does not unwind 
the soliton: the condensate remains in the soliton state with a phase jump of n. The phase jump or its unwinding 
could be detected by interference between two condensates, one of them in a ground state and the other with a soliton 
|p9| . A simpler way to verify that observed greying is due to the measurement induced unwinding is to change the 
measurement strength in a certain range and see if the greying time depends on the imaging laser intensity. We 
simulated the soliton unwinding for a range of measurement strengths. Fig. 5 shows the unwinding time r versus the 
laser intensity I. The time r is defined as that for which the density at x = in Fig. 4 achieves 10% of the maximal 
density. 




le-06 le-05 0.0001 0.001 0.01 0.1 

I ( mW / cm 2 ) 

FIG. 5. The unwinding time r as a function of the laser intensity 7. 



VII. CONCLUSIONS 



Quantum measurement theory has been applied to the dispersive imaging of a Bose-Einstein condensate. In 
the strong measurement limit the condensate is irreversibly driven into non-classical states with reduced number 
fluctuations. In the opposite limit of weak measurement the condensate can be approximately described for short 
timescales through a stochastic counterpart of the Gross-Pitaevskii equation. The latter has been applied for the study 
of the dynamics induced by dispersive imaging on a condensate prepared in a soliton state. The proposed model, 
besides allowing to intentionally design selective manipulation of the condensate state for instance to quench vortices 
without introducin appreciable depletion, could also lead to a better understanding of quantum phase transitions 
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in Bose condensates [h_2| . 
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